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Integrated  gasification  fuel  cell  (IGFC)  systems  that  combine  coal  gasification  and  solid  oxide  fuel  cells 
(SOFC)  are  promising  for  highly  efficient  and  environmentally  sensitive  utilization  of  coal  for  power 
production.  Most  IGFC  system  analysis  efforts  performed  to-date  have  employed  non-dimensional  SOFC 
models,  which  predict  SOFC  performance  based  upon  global  mass  and  energy  balances  that  do  not  resolve 
important  intrinsic  constraints  of  SOFC  operation,  such  as  the  limits  of  internal  temperatures  and  species 
concentrations.  In  this  work,  a  detailed  dimensional  planar  SOFC  model  is  applied  in  IGFC  system  analysis 
to  investigate  these  constraints  and  their  implications  and  effects  on  the  system  performance.  The  analysis 
results  further  confirm  the  need  for  employing  a  dimensional  SOFC  model  in  IGFC  system  design.  To 
maintain  the  SOFC  internal  temperature  within  a  safe  operating  range,  the  required  cooling  air  flow 
rate  is  much  larger  than  that  predicted  by  the  non-dimensional  SOFC  model,  which  results  in  a  larger 
air  compressor  design  and  operating  power  that  significantly  reduces  the  system  efficiency.  Options 
to  mitigate  the  challenges  introduced  by  considering  the  intrinsic  constraints  of  SOFC  operation  in  the 
analyses  and  improve  IGFC  design  and  operation  have  also  been  investigated.  Novel  design  concepts  that 
include  staged  SOFC  stacks  and  cascading  air  flow  can  achieve  a  system  efficiency  that  is  close  to  that  of 
the  baseline  analyses,  which  did  not  consider  the  intrinsic  SOFC  limitations. 

©  201 1  Elsevier  B.V.  All  rights  reserved. 


1.  Introduction 

Coal-based  power  plants  with  low  criteria  pollutant  emissions 
and  carbon  capture  capability  are  essential  for  satisfying  the  ever 
increasing  global  energy  demand  while  protecting  the  earth  from 
pollution  and  climate  change.  Conceptual  design  and  analyses  of 
the  integration  of  coal  gasifiers  with  high  temperature  fuel  cells 
have  attracted  much  research  interest  around  the  world  since  the 
early  1990s  and  the  results  have  clearly  shown  the  potential  of 
such  systems  for  superior  efficiency  and  emissions  compared  to 
other  system  approaches,  including  integrated  gasification  com¬ 
bined  cycle  (IGCC)  systems  [1-10].  Recently,  IGFC  systems  that 
take  advantage  of  more  efficient  catalytic  gasifiers  and  the  cooling 
caused  by  direct  internal  reforming  (DIR)  of  CH4  to  reduce  parasitic 
air  compression  load,  have  been  proposed  [11,12]  and  the  electri¬ 
cal  efficiency  of  such  systems  have  been  estimated  to  exceed  or 
approach  the  U.S.  Department  of  Energy  (DOE)  performance  goal 
for  next  generation  coal-based  power  plants:  producing  electricity 
at  60%  efficiency  (coal  HHV  basis)  while  capturing  more  than  90% 
of  the  inlet  fuel  carbon  in  a  pure  C02  stream  [13]. 
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However,  the  system  performance  of  these  designs  needs  fur¬ 
ther  scrutinization  because  most  of  the  published  efforts  are  based 
upon  use  of  non-dimensional  SOFC  models  that  only  analyze 
the  global  mass  and  energy  balances  of  the  SOFC.  As  an  elec¬ 
trochemical  energy  conversion  device,  the  performance  and  safe 
operation  of  SOFC  are  also  significantly  affected  by  the  internal  dis¬ 
tributions  of  temperature  and  species  concentrations,  which  are 
beyond  the  resolution  capability  of  non-dimensional  SOFC  mod¬ 
els.  To  clarify  this  issue,  a  detailed  dimensional  model  for  planar 
SOFC  has  been  previously  developed  for  use  in  such  IGFC  system 
analysis  work  [14].  The  current  work  demonstrates  and  discusses 
the  application  of  this  model  to  the  design  and  analysis  of  IGFC 
systems. 

To  provide  a  system  level  background  to  carry  on  the  discus¬ 
sions,  an  IGFC  system  consisting  of  catalytic  hydro-gasification 
and  hybrid  pressurized  SOFC  -  gas  turbine  power  block,  which 
is  one  of  the  most  promising  system  configurations  obtained  in 
previous  non-dimensional  SOFC  model  based  design  and  anal¬ 
ysis  [12],  has  been  chosen  as  a  “baseline”  case  in  this  work. 
The  detailed  dimensional  SOFC  model  is  introduced  in  the  place 
of  the  previously  used  non-dimensional  SOFC  model  and  the 
system  design  criteria  and  configurations  are  modified  accord¬ 
ingly  to  accommodate  the  new  constraints.  Differences  in  system 
design  approaches  and  system  performance  are  investigated  and 
presented. 
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Nomenclature 

A  area,  m2 

D/  eff  effective  diffusivity  of  species  i  in  porous  materials, 

m2  s-1 

E°  ideal  potential  of  H2  oxidization  at  ambient  pressure 

Fact  activation  energy,  kj  mol-1 

I<  convective  heat  transfer  coefficient,  W  m-2  K-1 

Kp  equilibrium  constant 

R  ohmic  resistance,  or:  heat  conduction  resistance, 

WK-1 

T  temperature,  K 

V  voltage,  V 

h  specific  enthalpy  of  species,  J  mol-1 

i  electric  current,  A 

j  electric  current  density,  A  m-2 

jo  exchange  current  density,  A  m-2 

p  pressure,  bar 

r  rate  of  reaction,  mol  s-1 

Xj  molar  fraction  of  species  i 

a  electron  transfer  coefficient 

y  pre-exponential  factor  in  exchange  current  density 

calculation 

8  thickness,  m 

Subscript 

IC  interconnect 

PEN  positive-electrolyte-negative  structure 

air  air  or  air  side 

amb  ambient  conditions 

an  anode 

cat  cathode 

ele  electrolyte  or  related  to  electrochemical  oxidation 

of  H2 

fuel  fuel  or  fuel  side 

rx  methane  reforming  reaction 

shift  water  gas  shift  reaction 

Superscript 
an  anode 

b  bulk  flow 

cat  cathode 

r  reaction  site 


2.  The  dimensional  SOFC  model 

Many  models  have  been  developed  to  provide  detailed  insights 
into  SOFC  operation  [15-18].  The  model  used  in  this  work  is  a 
quasi-two-dimensional  finite  volume  model  for  planar  SOFC  geom¬ 
etry.  The  model  has  been  developed  specifically  to  bridge  the  gap 
between  detailed  dimensional  SOFC  modeling  and  bulk  SOFC  mod¬ 
els  for  use  in  whole-plant  level  IGFC  system  analysis.  The  current 
model  thus  has  the  following  characteristics:  (1)  it  is  based  upon 
detailed  electrochemical  analysis  and  internal  heat  transfer  cal¬ 
culations;  (2)  it  can  give  not  only  SOFC  overall  performance  but 
also  internal  profiles  of  temperature,  current  density,  flow  com¬ 
positions,  etc.,  so  that  detailed  SOFC  operating  conditions  under 
different  system  configurations  can  be  investigated;  (3)  it  is  com¬ 
putationally  efficient  and  robust;  and  (4)  it  contains  the  flexibility 
to  be  directly  linked  to  system  analysis  tools  (such  as  chemical  flow 
sheet  software  Aspen  Plus®  [19]).  Special  attention  has  been  paid 
to  making  the  model  capable  of  reflecting  some  recent  develop¬ 
ments  in  SOFC  technology,  such  as  direct  internal  reforming  (DIR), 
anode  supported  geometry,  and  the  use  of  metallic  interconnects.  In 


Table  1 

SOFC  cell  single  channel  geometric  parameters  (co-flow  or  counter-flow). 


Flow  channel  length 

300  mm 

Flow  channel  width  (both  fuel  and  air) 

3  mm 

Fuel  channel  height 

1  mm 

Air  channel  height 

2  mm 

Anode  thickness 

1  mm 

Cathode  thickness 

0.05  mm 

Electrolyte  thickness 

0.01  mm 

Interconnect  thickness 

3.5  mm 

Rib  width 

2.42  mm 

addition,  the  electrochemical  parameters  of  the  model  were  deter¬ 
mined  through  sensitivity  analyses  to  match  the  SOFC  performance 
observed  by  modern  SOFC  developers  [14]. 

The  model  is  capable  of  simulating  the  two  parallel-flow  con¬ 
figurations  (co-flow  and  counter-flow),  which  are  believed  to  be 
sufficiently  representative  for  the  purposes  of  system  analysis. 
The  structures  of  fuel  flow  channel,  air  flow  channel,  positive 
electrode-electrolyte-negative  electrode  (PEN)  structure,  and  air- 
and  fuel-side  interconnects  (including  rib  structures)  are  resolved. 
The  geometric  configuration  of  the  model  is  shown  in  Fig.  1. 
The  anode-supported  configuration  that  minimizes  ohmic  losses 
through  use  of  a  very  thin  electrolyte  is  commonly  used  in  modern 
planar  SOFC.  Recently,  SOFC  producers  have  also  reported  success¬ 
ful  fabrication  of  cells  of  large  area  up  to  1000  cm2  (33  cm  x  33  cm) 
[20].  Throughout  this  work,  a  set  of  geometric  parameters  that  can 
well  represent  such  a  trend  of  planar  SOFC  size  and  performance 
has  been  employed  and  the  parameters  that  define  this  geometry 
are  listed  in  Table  1. 

The  fuel  cell  model  consists  of  two  interacting  modules:  the 
“species  conservation”  (SC)  module  and  the  “energy  conservation” 
(EC)  module.  The  SC  module  accounts  for  detailed  electrochem¬ 
ical  and  chemical  reactions  to  calculate  the  species  profiles  and 
current  density  distribution  within  the  SOFC.  It  is  assumed  that 
only  H2  participates  in  electrochemical  reactions,  while  CO  is  oxi¬ 
dized  through  the  water  gas  shift  reaction.  CH4  is  converted  into 
CO  and  H2  through  the  internal  steam  reformation  reaction,  which 
is  kinetically  limited.  Side  reactions  such  as  direct  electrochemi¬ 
cal  oxidation  of  CO,  or  dry  reforming  of  CH4,  are  assumed  to  be  of 
little  effect  and  are  not  included  in  the  calculations.  The  EC  mod¬ 
ule  calculates  temperature  distribution,  heat  transfer,  and  heat  loss 
throughout  the  fuel  cell.  Data  are  transferred  between  the  two 
modules  iteratively  until  the  calculation  convergence  criteria  are 
satisfied.  Some  important  equations  and  parameters  employed  in 
the  SC  module  are  summarized  in  Table  2.  Some  energy  conserva¬ 
tions  (for  co-flow  case)  employed  in  the  EC  module  are  summarized 
in  Table  3 ;  by  linearizing  the  species’  specific  enthalpies,  the  energy 
equations  in  the  EC  module  are  written  into  tridiagonal  matrices, 
which  can  be  solved  very  efficiently  by  the  tridiagonal  matrix  algo¬ 
rithm  (TDMA)  [21]. 

More  detailed  features  of  the  model  and  typical  modeling  results 
when  running  the  model  as  a  stand-alone  analysis  tool  were 
presented  in  a  previous  paper  [14].  In  addition,  the  model  has 
demonstrated  that  except  for  the  very  simple  case  of  hydrogen 
fuel  with  co-flow  configuration,  SOFC  operation  generally  results 
in  complicated  internal  temperature,  species,  and  current  density 
profiles  that  should  be  resolved  in  systems  analyses. 

This  work  focuses  upon  the  additional  insights  provided  by  this 
dimensional  model  (compared  to  a  non-dimensional  SOFC  model) 
in  an  IGFC  system  analysis  work  and  the  changes  in  system  design 
that  must  be  considered  due  to  the  insights  produced  by  the  dimen¬ 
sional  model.  To  achieve  this  goal,  the  developed  SOFC  model  was 
programmed  in  FORTRAN  language  and  linked  with  the  Aspen 
Plus®  process  engineering  flow  sheet  software  package.  A  user- 
defined  linking  interface  (also  programmed  in  FORTRAN)  is  set  up 
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Fig.  1.  Fuel  cell  geometry  for  co-  and  counter-flow  planar  SOFC  (cited  from  [14]). 


Table  2 

Summary  of  major  equations  and  parameters  employed  in  the  SOFC  species  conservation  module. 


Electrochemical  sub-model:  Vcen  =  VNernst  -  ^act  -  Oohm  -  Odif  =fU) 


Nernst  potential  [22] 

Activation  polarization  [23] 

Anode  side  exchange  current  density  jo, an  [14,24,25] a 
Cathode  side  exchange  current  density,  j0iCat  [14,24,25] a 
Ohmic  polarization 


Diffusion  polarization  [26] 


Water  gas  shift  reaction:  CO  +  H20  =  H2  +  C02 
Assumed  to  be  always  in  local  equilibrium  [25] 
Methane  reformation:  CH4  +  H20  -*  3H2  +  CO 

Kinetics  controlled  [27] 


VNernst=£0+^ 


VVi  Uo.51n(£L) 


where:  E°  =  1.28628053  -  2.8873  x!0-4TPE 


»ac.=  ^sinh-’(X) 

=  (B)  (£^)exP(-sSi)  where:  y3n  =  5.5  x  108  A  m-2  ;  £act,an  =  50  kj  mol-1 

Jo.cat  =  ycst ( B; ) ° 25 exp  (- j^r^j  where:  ya t-7  x  108  Am"2  ;  £act,cat  =  100  1<J  mol"1 


^7  ohm  “  l(^PEN  +  ^IC.fuel  +  ^IC,air) 


vh2oah2 


,  R-u  FpEN 

H - ln 

4F 


where  :  x(,  =  x\ 


jRu  FpEN^ai 


/VH20  -  aH20  2 FpanDar 

xo2=1+(xo2-1)exP 


h2  2  Fpa„Da„,eff 
jFijTpEN^an 


f  jRu  TpEN  ^cat  \ 
\  4FpcatDcat,eff  ) 


Ph2Pco2  _  xh2xco2 
Ph2oPco  —  xh2oxco 


=  exP  e-  3-961 


rrx  =  KrxPcH4  exp  J  A" 

where  :  yrx  —  4274  mol  s_1  m-2  bar-1 
Eact.rx  =  82  1<J  mol-1 


a  Parameters  obtained  through  sensitivity  analysis,  see  [14]  for  more  details. 


to  transfer  information  between  the  SOFC  model  and  the  Aspen 
Plus®  main  program  so  that  Aspen  Plus®  can  call  the  dimensional 
SOFC  model  automatically  in  an  iterative  calculation  process  and 
thus  eliminate  the  tedious  labor  and  possible  mistakes  that  might 
otherwise  be  involved  in  copying  and  transferring  large  amounts  of 
data  between  different  computer  programs.  Fig.  2  shows  the  repre¬ 
sentation  of  a  quasi-two-dimensional  co-flow  SOFC  model  in  Aspen 
Plus®,  where  the  streams  “FUEL”,  “OXID”,  “ANEX”,  and  “CAEX”  rep¬ 
resent  fuel  inlet,  oxidant  inlet,  anode  exhaust  and  cathode  exhaust, 
respectively;  “FCWORK”  is  the  electricity  produced  by  the  SOFC 


FUEL 

OXID 


COFLOW-1 


SOFC^ 


—n  ' 

I  1 


ANEX 

CAEX 


QLOSS 


FCWORI^ 


d- 

<■ 


Fig.  2.  Representation  of  a  quasi-two-dimensional  co-flow  SOFC  model  in  Aspen 
Plus®. 
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Table  3 

Summary  of  major  equations  employed  in  the  SOFC  energy  conservation  module,  co-flow  case. 


The  ith  control  volume  in  the  fuel  channel 


The  ith  control  volume  in  the  air  channel 


The  ith  control  volume  in  the  PEN 


The  ith  control  volume  in  the  interconnect 


y^nk(i)hk(i-  1)-  i  +  l)hfc(i)  +  i^fuel-^fuel  -  PEN  ( fpEN  ( 0  -  Tfuel(O)  + 

*WWi-ic(7jc(i)  -  Tfuei(i))  -  rrx(i)hcH4(i)  -  rrx(i)hH2o(0  +  rrx(i)h'co(i)  +  3 rrx(i)h'H2(i)  - 

rde(i)/iH2(i)  +  rele(i)h'H20(i )  =  0 

where:  k  =  H2,  CH4,  CO,  C02,  H20,  N2  and  Ar 

^  ^nk(i)hk(i  -  1 )  -  ^  ^Hk(i  +  1  )hfc(i)  +  /GirAair-PEN(7pEN(i)  -  Tair(* ))  +  fGirAair-Ic(Tc(l)  - 


k  k 

Tair(i))-ireie(i)ho2(i)  =  0 
where:  /<  =  02,  N2,  C02,  H20  and  Ar 


'rPEN(i-1)-'rPEN(i)  ,  rpEN(i+1)-rpEN(i)  ,  hctO-rpENCO 

rpen  rpen-ic 


rpen 

f^air-^air- PEN ( Fair (0  ~  T’pEn(i))  -  WeI e(0  +  rrx(l)hcH4 (i)  +  ^rx ( i )/1h2 0 ( 0 
3rrx(i)h'H2(i)  +  reie(i)hH2(i)  +  Jreie(i)ho2(i)  -  reie(i)h'H20(i)  =  0 


+  /'CfuelAfuel-PEN(’f’fuel(0  ~  Tpen(O)  + 

Wco(i)- 


MizlhMi)  +  ric(»+i)-hc(0  +  W)  +  Wfud-IC(  WO  -  Tic(i))  +  /CairAair_IC(Tair(i)  - 
T,c(i))  =  0 


stack  and  “QLOSS”  is  the  heat  loss  of  the  SOFC  stack  to  the  environ¬ 
ment.  SOFC  model  results  that  are  essential  for  the  system  analysis 
to  carry  on  (such  as  flow  rates  and  thermodynamic  properties  of 
SOFC  outlet  streams)  are  transferred  back  to  the  Aspen  Plus®  pro¬ 
gram  automatically,  the  rest  of  the  results  can  either  be  returned  to 
the  Aspen  Plus®  program  through  the  “User  Arrays”  form  or  saved 
in  a  separate  formatted  data  output  file. 

3.  The  “baseline”  case  and  performance 

The  “baseline”  IGFC  system  employed  in  this  work  is  the  most 
promising  design  obtained  in  previous  non-dimensional  SOFC 
model  based  IGFC  system  design  work  [12].  The  block  flow  diagram 
of  the  baseline  system  is  shown  in  Fig.  3.  The  system  consists  of  a 
catalytic  hydro-gasification  sub-system,  proven  low-temperature 
gas  cleaning  equipment,  and  a  hybrid  fuel  cell  -  gas  turbine 
power  block  (with  the  SOFC  operating  at  about  lObar).  The  sys¬ 
tem  uniquely  features  recycling  of  the  de-carbonized,  humidified 
SOFC  anode  exhaust  gas  back  to  the  hydro-gasifier  for  improved 
energy  integration. 

The  catalytic  hydro-gasifier  [28]  is  capable  of  producing  a  syngas 
with  high  CH4  content  at  high  cold  gas  efficiency.  The  CH4  con¬ 
tent  in  the  syngas  can  be  internally  reformed  in  the  SOFC  and  the 
chemical  cooling  provided  by  this  endothermic  reaction  can  signif¬ 
icantly  reduce  the  flow  rate  of  cooling  air  required  by  the  SOFC,  thus 
decreasing  the  parasitic  load  of  SOFC  stack  air  compression.  This  is 
a  very  important  contributing  element  to  the  overall  system  effi¬ 
ciency  of  several  recently  proposed  highly  efficient  IGFC  systems 
[11,12]. 

The  system  is  estimated  to  be  able  to  achieve  a  thermal  efficiency 
of  58.4%  while  capturing  94%  of  the  carbon  present  in  syngas  as  a 
compressed  CO2  stream.  The  efficiency  achieved  by  this  configu¬ 
ration  can  be  as  high  as  61.5%  if  it  is  not  required  to  compress  the 
separated  CO2  stream  to  a  high  pressure.  Detailed  system  perfor¬ 
mance  data  of  the  “baseline”  case  are  listed  in  Table  4. 

The  non-dimensional  SOFC  model  in  the  “baseline”  case  is 
replaced  with  the  dimensional  SOFC  and  the  system  performance 
is  scrutinized  in  detail,  as  will  be  shown  in  the  following  sections. 

4.  IGFC  system  analysis  using  dimensional  SOFC  models 

4.1.  Non-dimensional  SOFC  model  analysis:  a  recapitulation 

In  the  “baseline”  system,  because  the  non-dimensional  SOFC 
model  cannot  resolve  the  internal  temperature  and  species  dis¬ 
tributions  within  SOFC  stacks,  the  following  SOFC  performance 
characteristics  have  been  employed  or  assumed: 


1.  It  is  assumed  that  the  temperatures  of  the  fuel  and  air  flows  are 
very  close  to  those  of  the  solid  structures  (PEN  and  the  inter¬ 
connect  structure)  everywhere  inside  the  SOFC  cell,  thus  the 
temperatures  of  the  fuel  and  air  flows  are  good  approximations 
of  the  temperatures  of  the  solid  structures. 

2.  Fuel  and  air  inlet  temperatures  are  both  923  K  (650  °C),  this  can 
be  achieved  by  providing  sufficient  preheating  to  the  inlet  fuel 
and  air  flows. 

3.  It  is  assumed  that  fuel  and  air  outlet  temperatures  are  both 
1123  K  (850  °C). 

4.  A  design  criterion  requires  that  SOFC  maximum  operation 
temperature  does  not  exceed  1123K  (850  °C).  With  the  3rd 
assumption  described  above,  it  is  implicitly  assumed  that  the 
maximum  temperature  within  SOFC  channel  occurs  at  the  fuel 
and  air  outlet  sides. 

5.  Another  design  criterion  requires  a  maximum  SOFC  AT 
(difference  between  highest  temperature  and  lowest  SOFC  tem¬ 
perature)  that  is  no  larger  than  200  K.  Combining  the  2nd  and 
3rd  assumptions  described  above,  it  can  be  seen  that  the  maxi- 


Table  4 

Summary  of  thermal  performance  of  the  “baseline”  IGFC  plant. 


Major  gross  power  output  (MW) 

SOFC  electrical  power 

247.77 

Cathode  exhaust  expander 

63.38 

Steam  turbine 

2.53 

Reactor/expander  topping  cycle 

9.34 

Total  gross  power  generated 

323.26 

Major  auxiliary  power  consumption  (MW) 

Coal  milling  and  handling 

0.60 

Coal  pump 

1.60 

ASU  air  compressor 

1.55 

ASU  02  compressor 

0.55 

SOFC  air  compressor 

66.91 

Anode  exhaust  compressor 

4.33 

Recycled  H2  compressor 

3.91 

Selexol™  unit 

0.98 

MDEA  unit 

0.62 

Cooling  tower  fan 

0.45 

Cooling  tower  pump 

0.85 

Condensate  recycle  pump 

0.29 

BFW  feed  pump 

0.16 

Transformer  losses 

0.70 

C02  compression 

11.57 

Miscellaneous  BOP  &  lighting 

1.00 

Total  internal  power  consumption  and  losses 

96.50 

Overall  performance 

Net  electric  power  (MW) 

226.76 

Carbon  capture  rate 

93.7% 

Overall  thermal  efficiency  (%  coal  HHV) 

58.4% 
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Fig.  3.  Block  flow  diagram  of  “baseline”  case:  IGFC  power  only  plant  with  pressurized  SOFC  and  anode  exhaust  recycled  to  the  gasifier. 


mum  AT  is  implicitly  assumed  to  be  the  difference  between  the 
flow  outlet  temperature  and  inlet  temperature. 

6.  The  SOFC  maximum  air  utilization  factor  is  assumed  to  be  58%; 
this  will  approximately  result  in  02  mole  fraction  of  10%  in  the 
cathode  exhaust. 

A  non-dimensional  SOFC  model  confirmed  that  these  perfor¬ 
mance  characteristics  are  thermodynamically  correct  and  do  not 
violate  any  mass  or  energy  balance  requirements.  Flowever,  with 
the  insights  provided  by  the  detailed  dimensional  SOFC  models, 
it  is  found  that  some  of  the  assumptions  listed  above  turn  out  to 
be  unrealistic  or  erroneous,  i.e.,  impossible  to  achieve  given  the 
geometry,  physics,  chemistry  and  electrochemistry  extant  in  an 
SOFC. 


fuel  inlet  side  because  the  air  flow  has  picked  up  enough  heat  in 
the  cathode  channels  before  reaching  the  air  outlet  edge  (fuel  inlet 
edge)  and  thus  can  provide  enough  heating  for  the  endothermic 
steam  methane  reformation  occurring  at  the  location.  It  should  be 
pointed  out  that  the  counter-flow  configuration  can  still  encounter 
temperature  dips  that  are  severe  enough  to  hamper  operation 
when  the  syngas  contains  very  large  amounts  of  CH4  (more  than 
35vol.%  CFI4)  and  the  air  flow  rate  is  pushed  to  a  relatively  low 
level. 

Another  point  regarding  SOFC  operation  on  CFI4-rich  syngas 
is  that  the  counter-flow  configuration  generally  has  appreciably 
higher  average  current  and  power  density  than  the  co-flow  con¬ 
figuration.  This  finding  is  because  counter-flow  SOFC  with  internal 
steam  methane  reforming  can  result  in  species  and  temperature 
distributions  in  which  the  area  of  high  temperatures  (correspond- 


4.2.  Insights  provided  by  a  dimensional  SOFC  model 

SOFC  with  co-flow  and  counter-flow  configurations  in  IGFC  sys¬ 
tems  have  been  investigated  using  the  dimensional  SOFC  model 
described  above. 

It  is  found  that  when  the  inlet  fuel  contains  a  significant  amount 
of  CH4,  co-flow  operation  becomes  very  challenging  due  to  a  severe 
temperature  dip  at  the  fuel  and  air  inlet  edge.  In  a  co-flow  configura¬ 
tion,  at  the  fuel  inlet  edge  the  steam  methane  reforming  is  so  strong 
that  the  cell  temperature  drops  dramatically  and  the  chemical  and 
electrochemical  reaction  kinetics  are  all  slowed  down  significantly. 
The  cell  ends  up  with  very  low  fuel  and  air  utilization  as  well  as  very 
poor  current  and  power  density.  Fig.  4  shows  a  temperature  dip  in 
a  co-flow  configuration  SOFC.  The  SOFC  inlet  fuel  contains  about 
21.4%  H2,  4.0%  CO,  10.5%  C02,  30.9%  H20  and  27.6%  CH4  (molar 
fractions,  balanced  by  N2  and  Ar);  this  fuel  composition  is  typical 
in  the  IGFC  cases  discussed  in  this  work. 

In  contrast,  the  counter-flow  configuration  has  some  advantages 
and  is  generally  less  susceptible  to  severe  temperature  dips  at  the 
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Fig.  4.  Temperature  dip  in  co-flow  SOFC  operating  on  CH4  rich  syngas. 
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Fig.  5.  “Correct”  non-dimensional  SOFC  model  operation  conditions  result  in 
extremely  high  internal  temperature  spike  in  counter-flow  SOFC. 

ing  to  smaller  polarization  losses)  and  the  area  of  high  fuel/oxidant 
chemical  potential  difference  (leading  to  higher  Nernst  potential) 
overlap  each  other  in  relatively  larger  portion  of  the  cell  (see  [12] 
for  more  details  and  discussions  on  this  topic). 

In  IGFC  designs  in  which  CH4-rich  syngas  is  produced  to  fully 
take  advantage  of  the  cooling  effect  of  internal  steam  methane 
reforming  (such  as  the  “baseline”  case  in  this  work),  the  CH4  con¬ 
tent  in  the  fuel  to  the  SOFC  can  be  25vol.%  or  even  higher.  Thus, 
the  counter-flow  configuration  is  a  preferable  choice  in  these  cases 
and  worth  more  detailed  investigation. 

The  challenge  associated  with  a  counter-flow  configuration, 
however,  is  the  temperature  spike  that  occurs  within  the  SOFC. 
For  example,  using  a  quasi-two-dimensional  counter-flow  SOFC 
model,  it  is  possible  to  find  a  working  case  with  inlet  and  outlet 
flow  conditions  that  closely  match  those  working  characteristics 
obtained  from  a  non-dimensional  SOFC  model.  Fig.  5  shows  the 
internal  temperature  distributions  of  this  case.  It  can  be  seen  that 
the  temperature  rise  from  fuel  inlet  to  outlet  and  air  inlet  to  out¬ 
let  are  both  approximately  200  K,  which  matches  that  predicted 
by  a  non-dimensional  SOFC  model;  but  the  internal  temperature 
peak  in  the  PEN  structure  reaches  almost  1 700  K,  which  is  obviously 
much  higher  than  the  safe  operating  temperature  of  a  state-of- 
the-art  SOFC.  Thus,  while  thermodynamically  correct  operating 
conditions  can  be  defined  by  a  non-dimensional  SOFC  model,  actual 
performance  can  be  unacceptable  due  to  cell  internal  temperature 
distributions,  which  can  only  be  revealed  by  a  dimensional  SOFC 
model. 

4.3.  Required  modification  of  the  “baseline"  system 

Based  on  the  observations  described  in  Section  4.2,  when  putting 
the  quasi-two-dimensional  counter-flow  SOFC  model  into  the 
“baseline”  case  for  a  whole  system  analysis,  the  following  design 
criteria  must  be  added: 

1.  In  the  non-dimensional  “baseline”  design  the  air  flow  rate  to 
the  SOFC  stack  is  varied  to  achieve  about  50%  air  utilization. 
Now  to  mitigate  an  unacceptably  high  temperature  spike  in  the 
SOFC,  the  air  flow  rate  is  varied  so  that  the  AT  (defined  as  differ¬ 
ence  between  highest  temperature  and  lowest  temperature)  in 
PEN  structure  is  within  200  K.  In  this  case  the  actual  air  utiliza¬ 
tion  is  calculated  in  the  SOFC  model  and  reported  by  the  SOFC 
model  in  the  “User  Arrays”  form  in  Aspen  Plus®.  The  temperature 
of  the  PEN  structure  is  explicitly  reported  by  the  quasi-two- 
dimensional  model;  thus,  it  is  no  longer  necessary  to  estimate 
the  temperature  of  PEN  based  upon  the  temperatures  of  the 


fuel  and  air  flows.  It  is  implicitly  assumed  that  the  tempera¬ 
ture  of  PEN  is  very  close  to  that  of  the  interconnect  structure; 
this  is  justified  by  previous  dimensional  model  results  reported 
in  [12]. 

The  SOFC  fuel  utilization  is  determined  by  the  H2  requirement 
of  the  hydro-gasifier,  which  remains  the  same  as  in  the  “base¬ 
line”  case.  But  because  the  SOFC  geometry  is  now  resolved,  the 
desired  fuel  utilization  must  be  achieved  by  varying  the  number 
of  SOFC  stacks,  number  of  cells  in  each  stack  as  well  as  number 
of  channels  in  each  cell. 

The  inlet  fuel  and  air  streams  are  still  designed  to  be  preheated 
to  923  K  (650  °C),  this  design  will  help  maintain  a  reasonable 
cell  operating  temperature  level.  Together  with  the  1st  design 
criterion  described  above,  this  design  also  makes  sure  that  the 
peak  temperature  in  the  SOFC  does  not  exceed  1 1 23 1<  (850  °C). 
4.  The  SOFC  stack  chamber  (containment  vessel)  temperature  is 
varied  so  that  the  heat  loss  of  the  SOFC  stacks  is  about  1%  of  the 
heat  generated  during  chemical  and  electrochemical  reactions. 

The  SOFC  stack  working  voltage  is  set  to  be  0.8  V,  which  is  con¬ 
sistent  with  state-of-the-art  SOFC  performance  and  is  the  same  as 
the  value  assumed  in  the  “baseline”  case,  for  comparison. 

The  quasi-two-dimensional  counter-flow  SOFC  model  was 
plugged  into  the  “baseline”  case  to  replace  the  original  non- 
dimensional  SOFC  model  and  a  converged  Aspen  Plus®  simulation 
was  obtained,  which  is  designated  the  “single  stage  counter-flow 
SOFC”  case. 

Before  jumping  to  the  discussion  of  overall  thermal  performance 
of  the  modified  case,  it  is  important  to  first  observe  the  design 
changes  that  had  to  be  made  as  a  result  of  including  the  physics 
included  in  the  dimensional  SOFC  model.  Most  of  the  original  IGFC 
system  components  and  operating  designs  of  the  “baseline”  case 
remain  unchanged;  however,  the  important  changes  concern  the 
methanation  reactor  in  the  reactor/expander  topping  cycle  and 
the  SOFC  downstream  pre-reformer:  these  two  components  were 
removed  to  deal  with  the  much  colder  SOFC  anode  exhaust. 

In  the  IGFC  “baseline”  design,  because  the  pressure  of  the  clean 
syngas  is  much  higher  than  the  operating  pressure  of  the  SOFC, 
it  is  preferable  to  put  the  syngas  through  a  syngas  expander  to 
produce  some  useful  work  before  sending  it  to  the  SOFC  stack.  A 
shift  reactor  and  a  methanation  reactor  can  be  placed  before  the 
expander  to  increase  the  syngas  temperature  (both  the  water-gas- 
shift  reaction  and  the  methanation  reaction  are  exothermic)  for 
higher  overall  efficiency.  This  system  component  is  called  a  “reac¬ 
tor/expander  topping  cycle”  [29]  and  has  been  demonstrated  to  be 
able  to  increase  system  thermal  efficiency  appreciably  in  IGCC  sys¬ 
tems.  Because  of  the  methanation  reaction,  the  fuel  coming  out  of 
the  methanation  reactor  contains  a  large  amount  of  CFI4  (as  much  as 
36  vol.%  CH4).  Under  some  operating  conditions,  the  heat  required 
to  internally  reform  so  much  CFI4  is  even  larger  than  the  overall  heat 
released  during  the  chemical  and  electrochemical  reactions  in  the 
SOFC.  Such  a  global  energy  imbalance  around  the  SOFC  stack  can 
be  well  captured  by  the  non-dimensional  SOFC  model  (under  simi¬ 
lar  conditions  a  dimensional  SOFC  model  will  end  up  with  a  serious 
temperature  dip,  which  also  indicates  an  impractical  operating  con¬ 
dition).  Thus,  for  stable  SOFC  operation  some  of  the  CFI4  needs  to 
be  pre-reformed  externally,  this  is  achieved  by  splitting  some  of 
the  syngas  to  an  external  pre-reformer  with  the  remainder  being 
sent  to  a  pre-heater.  The  heat  requirements  of  the  pre-heater  and 
the  pre-reformer  are  both  provided  by  the  hot  SOFC  anode  exhaust. 
The  split  ratio  actually  determines  how  much  CH4  is  reserved  for 
internal  steam  methane  reforming  and  as  heat  sink  within  the 
SOFC  stacks.  In  the  “baseline”  case  a  design  specification  is  created 
in  Aspen  Plus®  to  automatically  vary  the  split  ratio  to  satisfy  the 
heat  balance  requirement  of  the  SOFC  stack.  The  reactor/expander 
cycle  design  is  schematically  shown  in  Fig.  6.  The  benefit  of  this 
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Fig.  6.  Reactor/expander  topping  cycle,  pre-reformer  and  pre-heater  design  in  the  “baseline”  case. 


design  is  that  the  sensible  heat  carried  by  the  anode  exhaust  can 
be  effectively  recovered  in  the  topping  reactor/expander  cycle  to 
produce  useful  work  and  the  system  becomes  more  efficient.  As 
described  above,  the  fuel  flow  outlet  temperature  is  assumed  to 
be  close  to  the  maximum  operating  temperature  of  the  SOFC  in 
the  non-dimensional  SOFC  model  analysis  and  thus  the  anode 
exhaust  is  at  high  temperature  which  enables  significant  heat 
recovery. 

On  the  other  hand,  when  the  quasi-two-dimensional  SOFC 
model  is  included  in  the  analysis  it  is  found  that  the  anode  exhaust 
temperature  is  not  as  high  as  that  predicted  previously.  The  internal 
temperature  profiles  of  the  SOFC  stack  in  this  new  working  case  are 
shown  in  Fig.  7.  It  can  be  seen  that  the  anode  exhaust  temperature  is 
about  943  K,  which  is  only  an  increase  of  20 1<  compared  to  the  fuel 
inlet  temperature  (instead  of  200  K).  This  small  temperature  rise 
would  make  it  very  challenging  for  heat  transfer,  requiring  very 
large  (and  expensive)  heat  exchangers  to  the  anode  inlet  stream. 
Moreover,  the  pre-reformer  is  more  energy  intensive  than  the  pre¬ 
heater  due  to  the  endothermic  steam  methane  reforming  reaction 
and  the  anode  exhaust  predicted  by  the  dimensional  SOFC  model 
no  longer  contains  sufficient  heat  energy  for  this  purpose;  thus,  the 


Fig.  7.  Internal  temperature  distributions  of  SOFC  in  the  “single  stage  counter-flow 
SOFC”  IGFC  design. 


pre-reformer  must  be  removed  in  the  new  design  that  employs  the 
dimensional  SOFC  model.  Correspondingly,  the  methanation  reac¬ 
tor  in  the  reactor/expander  topping  cycle  also  must  be  removed. 
Thus,  some  benefits  associated  with  the  reactor/expander  topping 
cycle  design  must  be  forfeited  when  considering  the  physics  of  a 
dimensional  SOFC  model. 

With  these  IGFC  system  differences  in  mind,  we  can  now  con¬ 
sider  and  discuss  the  thermal  performance  of  the  “single  stage 
counter-flow  SOFC”  case.  The  major  performance  characteristics 
of  this  new  case  are  compared  with  those  of  the  “baseline”  case  in 
Table  5,  along  with  the  performance  of  other  system  configurations 
that  will  be  discussed  later  in  this  paper.  It  can  be  seen  that  the  sys¬ 
tem  can  only  achieve  a  thermal  efficiency  of  42.1%,  which  is  more 
than  15  percentage  points  lower  than  the  result  obtained  from  the 
non-dimensional  SOFC  model  “baseline”  case.  This  is  mainly  due  to 
the  huge  parasitic  load  of  the  SOFC  air  compressor  that  is  required 
to  maintain  the  SOFC  within  the  temperature  constraints  (not  pri¬ 
marily  due  to  system  configuration  differences  discussed  above).  In 
Fig.  7  it  can  be  seen  that  the  SOFC  internal  peak  temperature  is  well 
controlled  to  be  lower  than  the  maximum  limit  and  this  is  achieved 
by  increasing  the  cooling  air  flow  rate  significantly:  the  non- 
dimensional  SOFC  model  predicted  that  the  SOFC  air  utilization  can 
be  as  high  as  0.58,  while  in  the  quasi-two-dimensional  SOFC  model 
analysis  the  air  flow  rate  has  to  be  two  and  a  half  times  larger  to  sat¬ 
isfy  the  cell  AT  and  peak  temperature  requirements,  resulting  in  a 
much  smaller  air  utilization  of  only  0.1 6.  It  can  be  seen  that  the  new 
system  has  a  much  larger  gross  power  output  thanks  to  the  larger 
air  flow  rate  to  the  cathode  exhaust  expander,  but  the  increase  in 
auxiliary  power  consumption  is  even  more  significant,  and  the  net 
effect  is  a  much  lower  net  power  output  and  lower  overall  system 
efficiency. 

Thus,  the  application  of  a  quasi-two-dimensional  SOFC  model 
has  revealed  more  stringent  requirements  regarding  the  opera¬ 
tion  of  SOFC  stacks.  To  satisfy  these  requirements,  particularly 
to  satisfy  the  peak  temperature  and  maximum  AT  require¬ 
ments  of  the  SOFC,  many  of  the  expected  benefits  associated 
with  internal  steam  methane  reforming  are  diminished,  result¬ 
ing  in  a  system  efficiency  that  is  significantly  lower  than  that 
previously  predicted  by  non-dimensional  SOFC  model  based 
analysis.  The  next  section  will  investigate  several  SOFC  mod¬ 
ule  design  options  that  may  be  used  to  still  achieve  high 
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Table  5 

Summary  and  comparison  of  the  thermal  performance  of  IGFC  systems  (the  “baseline”  case  uses  the  non-dimensional  SOFC  model,  while  all  other  cases  use  the  dimensional 
SOFC  model). 


Baseline 

Single  stage 
counterflow  SOFC 

Single  stage 
counterflow  SOFC 
(higher  voltage) 

Four  stage  air 
cascading 
counterflow  SOFCs 

SOFC  performance  parameters 

SOFC  operation  pressure  (bar) 

10.1 

10.1 

10.1 

10.1 

SOFC  voltage  (V) 

0.8 

0.8 

0.85 

0.8 

SOFC  single  pass  fuel  utilization 

0.73 

0.72 

0.73 

0.73 

SOFC  single  pass  air  utilization 

0.58 

0.16 

0.19 

0.46 

Fuel  outlet  temperature  (I<) 

1123 

943 

940 

970 

Air  outlet  temperature  (I<) 

1123 

994 

986 

1032 

Major  gross  power  output  (>1  MW) 

SOFC  electrical  power 

247.77 

247.31 

262.80 

247.83 

Cathode  exhaust  expander 

63.38 

178.60 

146.08 

72.10 

Steam  turbine 

2.53 

1.90 

1.90 

2.67 

Reactor/expander  topping  cycle 

9.34 

7.58 

7.58 

7.21 

Total  gross  power  generated  (MW) 

323.26 

435.62 

418.59 

330.41 

Major  auxiliary  power  consumption  (>1  MW) 

Coal  pump 

1.60 

1.60 

1.60 

1.60 

ASU  air  compressor 

1.55 

1.55 

1.55 

1.55 

SOFC  air  compressor 

66.91 

242.50 

203.33 

84.75 

Anode  exhaust  compressor 

4.33 

4.34 

4.34 

5.88 

Recycled  H2  compressor 

3.91 

3.94 

3.93 

3.91 

C02  compression 

11.57 

11.63 

11.63 

11.63 

Total  internal  power  consumption  and  losses  (MW) 

96.50 

272.14 

232.92 

115.96 

Net  electric  power  (MW) 

226.76 

163.49 

185.67 

214.45 

Overall  thermal  efficiency  (%  coal  HHV)  58.4% 

42.1% 

47.8% 

55.2% 

IGFC  system  efficiency  while  meeting  all  of  the  SOFC  design 
constraints. 


5.  Options  for  improving  IGFC  system  performance, 
considering  SOFC  species  and  temperature  distributions 

5. 1 .  Higher  SOFC  working  voltage 

A  higher  SOFC  working  voltage  is  desirable  because:  (1)  higher 
working  voltage  makes  the  SOFC  stack  more  efficient;  (2)  a  larger 
portion  of  the  chemical  energy  in  the  fuel  is  converted  into  electric 
power  and  less  is  released  as  heat.  Thus,  increasing  SOFC  operat¬ 
ing  voltage  may  help  mitigate  the  extreme  temperature  spike  and 
reduce  the  required  cooling  air  flow  rate  and  the  associated  com¬ 
pression  power.  While  this  design  change  could  benefit  and  could 
have  been  anticipated  with  a  non-dimensional  model,  the  signifi¬ 
cant  need  for  this  change  could  not  have  been  revealed  without  the 
dimensional  SOFC  model. 


With  the  progress  of  SOFC  performance  reported  in  the  recent 
years  and  the  expectation  that  SOFC  performance  will  continue  to 
improve,  a  high  working  voltage  of  0.87  V  has  been  assumed  in  the 
DOE  IGFC  report  analysis  work  [11].  In  the  “baseline”  IGFC  case,  a 
relatively  conservative  working  voltage  of  0.8  V  has  been  employed 
because  with  a  non-dimensional  SOFC  model  it  is  impossible  to 
determine  whether  or  not  the  SOFC  can  maintain  a  high  work¬ 
ing  voltage  due  to  either  increased  polarization  losses  or  reduced 
Nernst  potential  at  any  point  within  the  cell.  With  the  help  of  a 
detailed  quasi-two-dimensional  SOFC  model,  it  is  now  possible  to 
investigate  the  performance  of  an  IGFC  system  with  higher  SOFC 
working  voltage  while  considering  all  of  these  constraints. 

The  “single  stage  counter-flow  SOFC”  IGFC  case  was  modified 
and  the  cell  working  voltage  was  increased  from  0.8  V  to  0.85  V, 
the  system  thermal  performance  is  reported  in  Table  5  along  with 
that  of  the  “baseline”  case  and  the  “single  stage  counter-flow  SOFC” 
case. 

The  system  thermal  efficiency  improves  appreciably  from  42.1% 
to  47.8%,  thanks  to:  (1 )  relatively  higher  power  output  of  the  SOFC 
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fuel  4 
650  °C 
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cathode  exhaust 
753  °C 


anode 
exhaust  1 
671  °C 


anode 
exhaust  2 
689  °C 


anode 
exhaust  3 
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anode 
exhaust  4 
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Fig.  8.  Sankey  diagram  of  the  SOFC  module  design  with  four  stages  of  cascading  SOFC  stacks. 
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stack,  and  (2)  relatively  smaller  parasitic  power  consumption  of  the 
SOFC  air  compressor. 

However,  to  satisfy  the  SOFC  internal  temperature  require¬ 
ments,  a  large  amount  of  cooling  air  is  still  needed  and  the  SOFC  air 
utilization  is  only  0.19;  this  is  why  the  system  thermal  efficiency  is 
still  significantly  lower  than  that  predicted  in  the  “baseline”  case. 
Besides,  a  drawback  is  that  when  SOFC  working  voltage  is  increased, 
the  anode  exhaust  temperature  becomes  lower;  this  makes  it  more 
challenging  to  pre-heat  the  fuel  stream  with  the  anode  exhaust. 

It  should  be  noted  that  increasing  SOFC  working  voltage  gener¬ 
ally  results  in  smaller  average  current  density  and  for  a  prescribed 
total  power  output  more  SOFC  stacks  (and  higher  capital  cost)  will 
be  needed.  Thus,  to  determine  an  optimized  cell  working  voltage 
one  also  needs  to  take  into  account  economic  factors,  but  this  is 
beyond  the  scope  of  the  current  work.  Also,  it  should  be  noted  that 
due  to  the  theoretical  limit  of  Nernst  potential,  there  is  not  much 
space  to  further  increase  the  SOFC  working  voltage.  Thus,  although 
increasing  cell  working  voltage  can  indeed  improve  IGFC  system 
performance,  raising  the  voltage  alone  cannot  achieve  IGFC  system 
efficiency  as  high  as  that  predicted  by  the  non-dimensional  SOFC 
model  “baseline”  case,  while  considering  the  constraints  addressed 
by  the  dimensional  SOFC  model. 

5.2.  Cascading  of  identical  SOFC  stacks:  a  viable  design  strategy 

Based  upon  the  analysis  work  described  above,  it  can  be  seen 
that  a  key  factor  for  achieving  high  thermal  efficiency  is  to  reduce 
the  overall  cooling  air  flow  rate  (increase  air  utilization)  while  sat¬ 
isfying  the  SOFC  stack  safe  operating  temperature  requirements. 

An  SOFC  module  design  strategy  is  developed  herein  that  is 
comprised  of  cascading  identical  SOFC  stacks  with  series  air  flow, 
parallel  fuel  flow,  and  intra-stack  introduction  of  fresh  air  to  pro¬ 
duce  roughly  identical  operating  conditions  for  each  stack  in  the 


module.  This  design  may  be  able  to  achieve  high  IGFC  system 
efficiency  because:  (1)  air  utilization  in  each  stack  can  be  held 
low  enough  to  effectively  control  the  internal  SOFC  temperature 
peak;  (2)  the  cathode  exhaust  coming  out  of  one  SOFC  stack  still 
contains  large  amount  of  02,  which  can  be  used  in  downstream 
stacks;  (3)  typically  for  the  counter-flow  configuration  the  cathode 
exhaust  temperature  is  not  significantly  higher  than  the  air  inlet 
temperature,  which  assures  that  the  cathode  exhaust  coming  from 
upstream  SOFC  can  still  effectively  work  as  a  cooling  flow  in  down¬ 
stream  SOFC  stacks.  Plus,  fresh  makeup  air  can  be  added  in  after 
each  stack  to  provide  more  02,  to  cool  the  air  flow  to  the  desired 
temperature,  and  to  provide  roughly  identical  operating  conditions 
for  each  stack.  This  design  feature  of  the  developed  SOFC  stack 
module  concept  is  important  to  the  cost  of  the  system,  which  will 
depend  upon  the  mass  manufacturing  of  identical  stack  modules. 

With  the  dimensional  SOFC  model  successfully  linked  with 
Aspen  Plus®,  it  is  relatively  easy  to  investigate  various  SOFC  module 
configurations  and  parameters.  A  working  design  to  be  presented  in 
the  following  section  is  comprised  of  four  stages  of  cascading  SOFC 
stacks.  The  SOFC  sub-system  design  is  shown  in  Fig.  8,  with  the  plot 
representing  a  Sankey  diagram  [30]  where  the  width  of  the  streams 
represents  the  mass  flow  rate  of  each  stream.  The  working  voltage 
of  each  of  the  four  staged  SOFC  stacks  is  set  to  0.8  V.  The  total  incom¬ 
ing  fuel  flow  is  pre-heated  to  923  K  (650  °C)  and  distributed  evenly 
to  the  four  SOFC  stacks.  The  air  flow  is  cascaded  in  series  through  the 
four  stages  of  SOFC  stacks,  the  main  air  flow  fed  to  the  first  stage 
SOFC  is  also  pre-heated  to  923  K  (650  °C).  Small  amounts  of  cold 
fresh  makeup  air  are  mixed  with  the  cathode  exhausts  from  stages  I, 
II  and  III,  to  cool  down  the  main  air  flow  slightly  as  well  as  to  provide 
more  02  for  electrochemical  reaction  in  downstream  stages.  There 
is  no  need  to  pre-heat  the  makeup  air  because  relatively  lower  tem¬ 
perature  is  desired  for  cooling  purposes;  thus  makeup  air  coming 
out  of  the  air  compressor  at  about  603  K  (330  °C)  can  be  directly 
used.  In  the  Aspen  Plus®  simulation,  a  design  specification  is  set  up 
to  automatically  vary  the  main  air  flow  rate  so  that  the  AT  in  SOFC 
stack  I  (the  first  stage  SOFC  stack)  does  not  exceed  200  K;  three 
more  design  specifications  are  used  to  automatically  vary  the  three 
makeup  air  flow  rates  so  that  the  AT  in  stacks  II,  III,  and  IV  are  also 
no  larger  than  200  K. 

The  internal  PEN  temperature  distributions  of  the  four  SOFC 
stages  are  plotted  in  Fig.  9,  it  can  be  confirmed  that  the  PEN  tem¬ 
peratures  of  all  four  stages  of  SOFC  stacks  are  well  controlled.  More 
detailed  performance  results  for  the  four  SOFC  stacks  are  summa¬ 
rized  in  Table  6. 

It  should  be  noted  that  the  design  parameters  (geometry,  cell 
properties,  etc.)  of  each  of  the  four  stacks  have  been  kept  the  same 
(i.e.,  all  SOFC  stacks  in  this  cascaded  stack  design  are  identical)  and 
the  resulting  thermal  profiles  within  the  four  stacks  are  quite  sim¬ 
ilar  to  each  other.  This  will  ensure  that  the  same  exact  SOFC  stacks 
can  be  used  interchangeably  amongst  the  four  stages  in  this  design 
and  manufacturing  of  a  repeat  SOFC  stack  unit  will  help  reduce  the 
cost  of  the  SOFC  sub-systems. 


Table  6 

Performance  of  the  cascading  SOFC  stacks. 


Stack  I  Stack  II  Stack  III  Stack  IV  Four  stages 


Working  voltage  (V) 

Single  pass  fuel  utilization 
Single  pass  air  utilization 
Fuel  inlet  temperature  (K) 

Air  inlet  temperature  (K) 

Fuel  outlet  temperature  (I<) 

Air  outlet  temperature  (K) 
Maximum  PEN  temperature  (I<) 
Minimum  PEN  temperature  (K) 
AT  in  PEN  (I<) 


0.8 

0.8 

0.69 

0.72 

0.15 

0.16 

923 

923 

923 

942 

944 

962 

986 

1005 

1103 

1113 

907 

916 

196 

197 

0.8 

0.8 

0.74 

0.75 

0.17 

0.18 

923 

923 

961 

979 

977 

992 

1015 

1026 

1126 

1132 

927 

935 

199 

197 

0.8 

0.73  (overall) 

0.46  (overall) 

923 

N.A. 

970  (anode  exhausts  mixture) 
1026  (last  stage) 

N.A. 

N.A. 

N.A. 
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The  thermal  performance  of  the  whole  IGFC  system  with  the 
cascading  SOFC  stacks  module  design  is  presented  in  the  last  col¬ 
umn  of  Table  5,  along  with  the  performance  of  the  three  other  cases 
for  comparison.  By  employing  the  cascading  SOFC  module  design, 
the  overall  air  utilization  of  the  SOFC  stacks  can  be  increased  to  as 
high  as  0.46;  this  is  much  closer  to  the  “ideal”  prediction  made  in 
the  “baseline”  case  with  a  non-dimensional  SOFC  module.  The  par¬ 
asitic  power  consumption  of  the  SOFC  air  compressor  also  drops 
significantly,  thanks  to  the  much  smaller  overall  cooling  air  flow 
rate.  The  system  manages  to  achieve  a  thermal  efficiency  of  55.2%, 
which  is  only  about  3  percentage  points  lower  than  that  predicted 
using  the  non-dimensional  model  in  the  baseline  case.  This  overall 
IGFC  system  efficiency  is  very  impressive  considering  all  the  strin¬ 
gent  temperature  control  requirements  that  have  been  imposed, 
and  is  especially  promising  in  comparison  to  other  strategies  for 
electricity  production  from  coal  [31  ].  Besides,  the  last  stage  cathode 
exhaust  temperature  and  the  anode  exhaust  mixture  temperature 
are  both  significantly  higher  than  those  achieved  in  the  single  stage 
SOFC  cases,  which  is  beneficial  to  the  IGFC  system  design  and  oper¬ 
ation  of  the  pre-heaters. 

A  down  side  of  the  design  is  that  some  pressure  head  of  the 
compressed  flows  is  wasted  because  of  the  cascading  stages,  which 
is  also  part  of  the  reason  why  the  system  overall  efficiency  is 
still  lower  than  that  predicted  in  the  non-dimensional  baseline 
case. 

Overall  the  performance  of  the  IGFC  with  cascading  SOFC  stacks 
is  considered  satisfactory  and  is  the  most  promising  system  iden¬ 
tified  thus  far. 

6.  Summary  and  conclusions 

This  work  investigates  the  application  of  a  detailed  dimensional 
planar  SOFC  model  in  the  design  and  analysis  of  coal-based  IGFC 
systems.  A  detailed  quasi-two-dimensional  finite  volume  SOFC 
model  [14]  has  been  linked  with  Aspen  Plus®  to  perform  the 
analyses.  The  most  promising  case  identified  in  previous  non- 
dimensional  SOFC  model  based  design  and  analysis  work  [12]  has 
been  used  as  a  starting  point  and  “baseline”  case. 

Two  major  differences  are  observed  when  the  dimensional 
SOFC  model  is  used  in  place  of  the  original  non-dimensional  SOFC 
model:  (1)  the  anode/cathode  exhaust  temperatures  end  up  much 
lower  than  those  predicted  using  the  non-dimensional  SOFC  model, 
which  makes  the  heat  transfer  in  the  pre-heaters  more  challenging 
and  renders  the  methanation  reactor  and  external  pre-reformer  in 
the  original  “baseline”  design  no  longer  useful;  (2)  the  cooling  air 
flow  rate  required  to  control  the  internal  temperature  spike  and  AT 
is  much  larger  than  that  predicted  by  the  non-dimensional  model, 
which  significantly  increases  the  parasitic  power  consumption  of 
SOFC  air  compressor  and  seriously  diminishes  the  system  thermal 
efficiency.  The  modified  system  that  only  considers  the  additional 
constraints  imposed  by  the  dimensional  SOFC  model  has  an  effi¬ 
ciency  of  42.1%,  which  is  disappointingly  low  for  the  novel  IGFC 
system  (even  though  it  is  still  much  higher  than  comparable  IGCC 
cycles). 

Design  changes  to  improve  the  overall  IGFC  system  efficiency 
while  fully  accounting  for  the  SOFC  operating  constraints  using 
the  dimensional  model  have  also  been  investigated.  By  increasing 
the  SOFC  working  voltage  from  0.8  V  to  0.85  V,  the  IGFC  system 
thermal  efficiency  increases  from  42.1%  to  47.8%,  this  is  an  impres¬ 
sive  improvement  but  the  system  performance  is  still  significantly 
lower  than  that  predicted  by  the  non-dimensional  SOFC  model.  A 
design  where  air  flow  is  cascaded  in  series  through  four  stages  of 
SOFC  stacks  is  found  to  be  quite  an  effective  solution.  In  this  design, 


multiple  stages  of  identical  SOFC  stacks  are  employed,  with  paral¬ 
lel  fuel  flow  and  series  air  flow  (with  some  fresh  makeup  air  added 
after  each  stage).  The  design  can  achieve  low  air  utilization  (suffi¬ 
cient  cooling)  within  each  stage  while  maintaining  relatively  high 
overall  air  utilization.  An  SOFC  sub-system  comprised  of  four-stage 
cascading  stacks  can  achieve  an  overall  IGFC  system  efficiency  of 
55.2%,  which  is  much  closer  to  the  “ideal”  performance  (58.4%) 
predicted  using  the  non-dimensional  SOFC  model. 

This  work  clearly  demonstrates  the  advantages  and  importance 
of  using  detailed  dimensional  SOFC  models  to  avoid  over-simplified 
or  even  erroneous  SOFC  operating  assumptions  in  IGFC  systems 
analyses.  Future  work  will  explore  the  application  of  the  dimen¬ 
sional  SOFC  model  in  other  IGFC  system  configurations,  such  as 
cases  with  H2  co-production  or  those  that  use  oxygen-blown  gasifi¬ 
cation,  in  which  detailed  and  accurate  knowledge  of  the  SOFC  stacks 
is  also  essential  to  the  design  work. 
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